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Abstract 

We study decidability and complexity questions related to a continuous analogue 
of the Skolem-Pisot problem concerning the zeros and nonnegativity of a linear 
recurrent sequence. In particular, we show that the continuous version of the 
nonnegativity problem is NP-hard in general and we show that the presence 
of a zero is decidable for several subcases, including instances of depth two or 
less, although the decidability in general is left open. The problems may also be 
stated as reachability problems related to real zeros of exponential polynomials 
or solutions to initial value problems of linear differential equations, which are 
interesting problems in their own right. 

Key words: Skolem-Pisot problem. Exponential polynomials, Continuous 
time dynamical system. Decidability, Ordinary differential equations 



1. Introduction 

Skolem's problem (also known in the literature as Pisot's problem) asks 
whether it is algorithmically decidable if a given linear recurrent sequence (LRS) 
has a zero or not. A LRS may be written in the form: 

Uk = On-lUk-l + a„_2Ufe-2 H h aoUk-n, 

for k > n where uo, ui, . . . , Un-i G Z are the initial inputs and oo, ai, . . . , On-i S 
Z are the recurrence coefficients, see also [9]. This forms the infinite sequence 
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(■"fc)fe°=o — ^- assume ao is nonzero, otherwise a shorter and equivalent 

recurrence exists. Such a recurrence sequence is said to be of depth n. 

For a hnear recurrent sequence u = {uk)'^Q C Z the zero set ofu is given by 
Z{u) = {i G N\ui = 0|. One of the first results concerning the zeros of LRS's 
was by T. Skolem in [23|, when he proved that the zero set is semilinear (i.e., 
the union of finitely many periodic sets and a finite set). This result was also 
later shown by K. Mahler [l8| and C. Lech 17 1 and is now often referred to as 
the Skolem- Mahler- Lech theorem. It is known that determining if Z[u) is an 
infinite set is decidable as was proven by Berstel and Mignotte Q . 

It was shown by N. Vereshchagin in 1985 that Skolem's problem (i.e., the 
problem "is the zero set of a LRS empty?" ) is decidable when the depth of the 
linear recurrent sequence is less than or equal to four in [2^ . It was also recently 
shown that Skolem's problem is decidable for depth five in [l3|, but the general 
decidability status is open. It is also known that determining if a given linear 
recurrent sequence has a zero is NP-hard, see [gj. 

Note that we may always encode a linear recurrent sequence of depth n into 
an integral matrix A G z("+i)^("+^) such that Uk = A\^^^-y for fc > 1. This 
follows since given the initial vector u — (uqi ui, • ■ • , Un~i)^ and the recurrence 
coefficients, ao, ai, . . . , a„_i, we first define matrix A' G Z"'*": 



A' 



Note that {A'^u = {uk,Uk^ 
1 dimension to give: 
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, Uk+n-i)- Now we shall extend this matrix by 



A'v 




where is the zero vector of appropriate size. It is not difhcult to now see that 
Uk — Ai ^ for fc > 1 as required. Skolem's problem in this context is therefore 
to determine if the upper right entry of a positive power of an integral matrix is 
zero. More generally, one can show that Skolem's problem is equivalent to the 
following problem: given a matrix A G F"^" and two vectors c, xq G F" , is there 
a nonnegative integer t such that c^A'^xq = 0? We add that a generalization 
of this problem where we may take any product of two integral matrices of 



dimension 10 is known to be undecidable, see 

In this paper we shall consider a dynamical system whose updating trajectory 
is given by = Ax{t) where A G W^" and the initial point a;(0) G M" is 
given. We shall be interested in determining whether this trajectory ever reaches 
a given hyperplane, thus the problem is equivalent to determining if there exists 
t G K>o such that exp{At)x{Q) — where c G M" defines the hyperplane. We 
consider this as the Skolem-Pisot problem in continuous time. We show that 
for instances of size two or less this problem is decidable. 
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We shall also show that determining if exp(At)a;(0) reaches zero is com- 
putationally equivalent to determining whether a given real-valued exponential 
polynomial f{z) = X^jLi ^'^P(^i-^)i where each Pj is a polynomial, ever 
reaches zero for a positive real value. This is also equivalent to determining if the 
solution y{t) of an ordinary differential equation y^'^^ + ak-iy^'^~^^ + ■ ■ ■ + aoy = 
with given initial conditions y'''^~^^ {0),y^'^~'^\0), . . . ,y{0) ever reaches zero. 

From 1920, Polya and others characterized the asymptotic distribution of 
complex zeros of exponential polynomials 0, [2l|, [g^, [2J, [2^ HI] • Upper bounds 
were also found on the number of zeros in a finite region of the complex plane, 
using the argument principle. Less is known about real zeros. Upper and lower 



bounds on the number of zeros in a real interval are given in |27| . A formula 
for the asymptotic densit y o f real zeros for a restricted class of exponential 
polynomials was found in [l5| . Some observations on the first sign change of a 
sum of cosines are collected in [i^l- However, no criterion has been proposed to 
check the existence of a real zero for a real exponential polynomial. 

A related problem, determining whether a given linear recurrent sequence 
has only nonnegative terms, the nonnegativity problem, is decidable for dimen- 



sion 2, see [1^. The authors note that if the nonnegativity problem is decidable 
in general, it implies Skolem's problem is decidable. This follows since if (ufc)^Q 
is recurrent, then so is {u\ — l)^o- 

We may note that using the linear recurrent sequence (ufc)^Q from the 
proof of NP-hardness of Skolem's problem in ^ , and converting it to the form 
(m^. — 1)^0, allows one to easily derive the following result: 

Theorem 1. It is NP-hard to decide if a given linear recurrent sequence is 
nonnegative, i.e., the nonnegativity problem is NP-hard. 

This holds since if (ufc)^o is represented by a matrix Z"^", then {uf. — l)^g 
may be represented by a matrix Z^" +i)x(n +i) thus we have a polynomial 
time reduction. In this paper we show that the nonnegativity problem in the 
continuous setting is also NP-hard. 

Given a matrix M G R"^" and vectors u,v £ R", the orbit problem asks if 
there exists a power /c S N such that M^u = v. Thus it is a type of reachability 
problem, see [7[. This was shown to be decidable even in polynomial time, 
see The corresponding version of this problem for continuous time asks 

whether for a given M G M"^" and vectors a,b € R" there exists some t € ]R>o 
such that exTp{Mt)a = b. This problem was proved to be decidable in 



2. Preliminaries 

Let A e F"^" denote an n x n matrix over the field F and cr(A) the set of 
eigenvalues of A. For a complex number z G C we denote by 5R(z) the real part 
of z and by 3(z) the imaginary part of z. We use the notation M>o to denote 
the nonnegative real numbers. 

We shall denote an exponential polynomial / : C — > C by a sum of the form: 
f{z) = T,T=i ^j (^) exp(6'jz), where Pj e C[X] and 9j e C. 
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Given a matrix A G C"^" we shall denote by the dominant eigenvalues of 
A the set of eigenvalues of A with maximum real part, i.e., 

{6 e <7{A)\^{e) > ^{9'),e' e 

We will later require the following theorem from Diophantine approximation 

i: 

Theorem 2. (^Baker^ Let ai, . . . , Offe, /Sq, . . . , be algebraic numbers. Then 
the combination 

A = /3o + ^ /?, In ai 

i 

is either zero or satisfies |A| > h^^ , where h is the largest height of (3i, . . . ,Pk, 
and N is a computable constant depending only on Inai, . . . , Ina^ and the max- 
imum degree o/ /?o , • ■ • , /3fc • 

Recall that for an algebraic number /3 with minimal polynomial 

0<i<!i 

its degree is d and its height is max \ai\. We shall also use the following theorem 
regarding the transcendence degree of the field extension of algebraic numbers 
when considering their exponentials: 

Theorem 3. (^Hermite-Lindemann^ - Let a^, Xj e C for < j < n — 1 be 
algebraic numbers such that no aj — and each Xj is distinct. Then: 

n~l 

j=o 

The following theorem concerns simultaneous Diophantine approximation of 
algebraic numbers which arc linearly independent over the rationals. 

Theorem 4. ('Kronecker, see [8^]) Let 1, Ai, A2, . . . , A„ £ M &e real algebraic 
numbers which are linearly independent over Q. Then for any a £ M" and 
e > 0, there exists p S Z" and /c e N such that \{kXi ~ ai ~ pi)\ < e for all 
1 < i < n. 

3. Skolem's Problem in Continuous Time 

We shall consider continuous time systems governed by the rule '^^^p = Ax{t) 
where ^ is a real matrix and x(t) is a real vector 0. We are interested in the 
decidability of whether from an initial vector a;(0), we cross a given hyperplane. 



^We consider entries to be algebraic so that the input to a problem has a finite description. 
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We may consider this as a "point-to-set" reachability problem in a dynamical 
system, see 0] for other examples. 

Let = Ax{t) where A e R"^" and x{t) e M". Given the initial vector 
a;(0) G M", then x{t) is given by: 



oo 



xit) = exp(At) • a;(0) ^ • a;(0). 

j=o -^^ 

Given a vector c e M" defining a hyperplane, we would like to determine 
if there exists some t e R>o such that c^x{t) — 0. In other words, whether 
the flow of the point x(0) ever intersects the hyperplane. If such a t exists, we 
say that there exists a solution to the instance q An instance of Continuous 
Skolem Problem therefore consists of the matrix ^ e R"^", the initial point 
x{0) S R" and the hyperplane vector c £ M". 

3.1. Equivalent Formulations 

To analyze the behaviour of the system, we will convert a given instance 
of Continuous Skolem Problem into various forms which have different 
properties but which are essentially equivalent to the original problem. 

Given such an instance, the following lemma shows that the problem is 
equivalent to determining if the upper right entry of the exponential of a matrix 
equals some constant real. A similar construction is known in the discrete case 
as shown in Section [TJ 

Theorem 5. Given an instance of CONTINUOUS Skolem Problem defined 
by f{t) = c^exp(Ai)a;(0) where A G R"""" and c,x(0) G R". There ex- 
ists a polynomial-time computable matrix B G R("+2)^("+2) such that f{t) = 
exp(_Bi)i „_|_2 + A, where A — c^a;(0) G R is constant. 

Proof. We are given the function f{t) = exp(At)x(0). Let B G m("+2)x("+2) 
be given by: 

/ c^A c^Aa;(0) 
5=0 A Ax{Q) 
V 0^ 

where = (0, 0, . . . , 0)"^ G R", thus: 

' exp(A) - exp(A)a;(0) - A 

exp(B) = ( exp(A) exp(A)x(0) - a;(0) 

0^ 1 



where A = C"'"x(0) is constant. This can be seen from the power series represen 

jj=0 j 



tation exp(ii3) = J27Lo Tf^'' ■ Therefore f{t) = exp{i?t)i_(„+2) + A and thus an 



^Note that, in the style of Skolem's problem, we shall be more interested in determining 
whether any solution exists, rather than trying to find an algebraic description of the solution. 
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instance of Continuous Skolem Problem can also be given by a single real 
matrix B and the problem of whether f{t) reaches zero for t £ M>o is equivalent 
to whether exp(St)i („_|_2) ever equals —A. □ 

Theorem 6. The following problems are computationally equivalent with poly- 
nomial time reductions (where all parameters are algebraic numbers): 

(i) Does there exist a solution to a given instance of CONTINUOUS Skolem 
Problem? 

(ii) Determine if a real-valued exponential polynomial: 

m 

f{t) = Y,PAt)e'>^\ 

has a nonnegative real zero (where 0j € C and Pj G C[X]j. 
(Hi) Determine if a function of the form: 

rn 

fit) = E ^"^''(P^At) cos(A,t) + P^At) sin(A,t)) 

has a nonnegative real zero (where rj,Xj £ R and Pi j G R[X]J. 

(iv) Determining whether the solution y{t) to an ordinary differential equa- 
tion y'^^^ + ak-iy'^^~^^ + . . . + a^y — with the given initial conditions 
y('=-i)(^0), y^'^~^^(0), . . . , y(0) reaches zero for a nonnegative real t. 

Proof, [i) => (ii): Let J e C"^" be the Jordan matrix for A, thus we may write 
A = PJP-^ for some P e GL{n,C)E Since exp(PJP-i) = Fexp(J)p-i, we 
can ask the equivalent problem, does there exist a time t > at which: 

c^y(t) = c^exp(tA)y(0) 
— exp{tJ)v — 0, 

where u, w G C" are defined by = c^P and v = P~^y{Q)l 

Let J = Ji© J2®. . .® Jm be a decomposition of J into a direct sum of Jordan 
blocks with Ji G C"iX"i and X^i^li = Each Jordan block may be written 
Ji = Oiljii + Mi where 9i E C is the associated eigenvalue, 1^ G Z"'^"' is the 
identity matrix and Mi G Z"'^"' has 1 on the super-diagonal and elsewhere. 

For 1 < i < m, we see that 9il„- and Mi commute and therefore exp(t J^) = 
exp(i6'i/„. ) exp(iMi). The value of exp(i^?i/„J is e*^'J„.. Let exp{tMi) = 
[rrijk] G Q"^", then 

f • if , < fc 

I ; otherwise 



•^These can be effectively found since we only need algebraic descriptions of the Jordan 
normal form J and the similarity matrix P. 
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Therefore we may convert our problem equivalently into deciding whether 
there exists a. t € R>o such that f{t) = where / : M ^ C is defined by: 

rn 

fit) = u^eMJt)v = 5]P,(t)e^^*, (2) 

and Pj G C[X] are polynomials and whose degree depends upon the size of 
the corresponding Jordan block and 9j G C. The polynomials Pj can be de- 
rived from Equation ([T|). Note that each of these steps is effective and can be 
computed in polynomial time for algebraic entries of the initial matrix A. 

(ii) [Hi): This results from Euler's formula for the complex exponential 
and the fact that f{t) is a real valued function. 

(Hi) (iv): Functions of the type 

rn 

f{t) = J2 ^"^''(Phjit) cos(A,i) + P2,jit) sin(A,i)) 

where rj,Xj £ M are fixed and Pkj are arbitrary real polynomials of degree 
< dj form a real vector space of dimension k = "^Y^^^^idj + 1). This vector 
space is closed under differentiation. Hence the first (fc + 1) derivatives of / 
are related by f'^'^^ + afe-i/^*^^^) + . . . + oq/ = where each aj can be found 
in polynomial time. By Cauchy's theorem for ordinary differential equations, 
a function / is completely determined by the given relation and the initial 
conditions /C^-^) (0), /C'-^) (0), . . . , /(O). 

(iv) => (z): The characteristic equation of the linear homogeneous differential 
equation is given by z'' + z^~^ak-i + . . . + oq = 0. It is well known that we can 
form the companion matrix of the equation in order to convert the problem into 
an instance of CONTINUOUS Skolem Problem. The initial values are then 
present in the initial vector a;(0). □ 

Lemma 7. Let A e M"^" and c, a;(0) £ M" form an instance of Continuous 
Skolem Problem. For any \ eC we may form a system fx{t) = exp{t{A+ 
XI))v where u,v ^ C", <j{A + XI) = (t{A) + X and f{t) ~ if and only if 
hit) = 0. 

Proof. Let A G C and define y{t) — e'^*x(t), thus: 

Ml - Xe^tx(t) + e^'^^ 
dt " ^'^ dt 

= e^\XI + A)x{t) 
= {XI + A)y{t) 

Define Ax = XI + A, thus: 

y{t) = exp{tAx)y{0). 

Note that there exists t > such that c^x{t) = if and only if c^y{t) = 
0. □ 
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As an example, which will be useful later, let us set A = —max{^{0)\0 G 
a{A)}, so that all eigenvalues are shifted to the left complex half-plane or the 
imaginary axis. This means that we have, in effect, split the set of eigenvalues 
into two sets, one which decays exponentially with time and one which consists 
of purely imaginary values. 

We now remark that any nontrivial solution to the problem will in fact be 
transcendental. 

Theorem 8. Given an instance of CONTINUOUS Skolem Problem, all so- 
lutions, if any exist, are transcendental unless the polynomials Pj{t) share a 
common positive real root. 

Proof. The corresponding exponential polynomial formed as in Theorem [6] will 
be in the form: 

m 

We may assume no Pj G C[X] is zero otherwise simply remove it from the sum 
and that each 9j is distinct, otherwise group them together. Thus, according 
to Theorem [3] (the Hermite-Lindemann theorem), this exponential polynomial 
only has solutions for transcendental times t where t G K>o . □ 

4. Decidable Cases 

We shall now investigate some classes of instances for which CONTINUOUS 
Skolem Problem is decidable. 

Theorem 9. The Continuous Skolem Problem for depth 2 is decidable. 

Proof. Assume we have an instance of Continuous Skolem Problem given 
by f{t) = (ci,c2)exp(At)(xi,X2)'^ with A e M?""^. Let S G GL(C,2) put A 
into Jordan canonical form. We can rewrite f{t) = (ai, (32) exp(Ji)(/3i, /?2)"^, 
where J = S~^AS is a Jordan matrix. 

If A has one eigenvalue 9, with algebraic multiplicity 2, then 6* G M. If 
9 has geometric multiplicity 1 then by Theorem [5] we must solve an equation 
of the form (1 + xt)ye*^ where x, 1/ G M, thus the instance has a solution if 
and only if — 5 G K>o. If 9 has geometric multiplicity 2 then we must solve 
e^*(ai/3i + a2/92) = which has a solution if and only if (ai/3i + 02/32) = 0. 

Otherwise, J is diagonal and we must determine if there exists a i G M>o 
such that e*''i + 06*"^ = for a G R. Either 6*1, 6I2 G K or = ^ G C. 

li 9i,92 G K. assume without loss of generality that 9i < 92 and we have 
g(t) = e*^i + ae*^^ thus, by taking logarithms, t — ^'^j^Zg.^ is a solution of 
g(t) = and thus there exists a solution if and only if ^^j^Zg^ ^ ^>o- 

In the other case 9i =92 G C. Since we may therefore shift the real part as 
allowed by [71 assume that 6'i, 6*2 G iM. At time t = 2^3(01 ) have 

= cos(|) + acos(-f ) = 
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which is a solution, thus we are done. 



□ 



The following theorem shows that the class of instances where all elements of 
the input are nonnegative reals in the continuous setting is trivially decidable in 
polynomial time, whereas in the discrete time case, the problem is NP-hard, as 
shown in Q. In fact, using Lemma [7l we see that in the continuous setting the 
Skolem-Pisot problem is polynomially decidable even where the matrix given is 
a Metzler matrix^ meaning only off-diagonal elements need be nonnegative. 

Theorem 10. For an instance of CONTINUOUS Skolem Problem given by 
A G M"^" and c, a:(0) G R"q where A is a Metzler matrix (thus all off-diagonal 
elements are nonnegative) and fit) = exp(j4t)a;(0), then we may decide if 
there exists a solution in polynomial time. 

Proof. Let A be the minimal diagonal element of ^. If A < then by Lemma [71 
we may form an equivalent instance A' ~ A + \I where A' e M"^". Thus 
assume without loss of generality that A is a nonnegative matrix and c, a;(0) are 
nonnegative vectors. 

Note that exp(i2A) > exp(tiA) for any ^2 > £ which is a consequence 
of the power series representation of exp(j4t) ~ Sjlo ^'^'^ ^^^^ thaX 
A e Rl"^''. We see that /(O) = c^x{0) e R>o. Now, if /(O) = then this 
is a solution, otherwise, since the matrix exponential increases monotonically 
componentwise with time for a nonnegative matrix, there exists no solution. □ 

In some special cases, some eigenvalues of A do not influence the function 
f{t). This is easily seen when A is put in its Jordan form J = P^^AP: 



where m, u G C" are defined by = c'^P and v — P^^x{0). Obviously, if the 
entries of c or x(0) corresponding to a particular Jordan block are zero, this block 
does not play any role and one may remove it without changing the function 
f{t). More generally, from Equation ([3]) it is easy, as shown in Theorem [6l to 
write the function / as follows: 



where the 9j are the distinct eigenvalues of A and the Pj are complex poly- 
nomials. If no polynomial Pj{t) is identically zero, then we say that the triple 
(A, c, x{0)) is reduced. If some of the Pj are zero, we can remove the correspond- 
ing terms from Equation ([4]) , since it does not change the value of /. 

Theorem[6] shows how to build an equivalent instance of the form {A' , c' , x' {Q)) 
from an instance of the form of Equation ([4]) . One would then obtain a reduced 
instance of the Continuous Skolem Problem. This can be done in a pre- 
processing phase. 



f{t) = u'^exp{tJ)v 



(3) 




(4) 
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Theorem 11. Let = Ax{t) for A e M"""" andx{t) £ M" define an instance 
of Continuous Skolem Problem given by f{t) = exp{At)x{0) = 0. // 
{A,c,x{0)) is reduced and none of the dominant eigenvalues of A are real then 
the problem is decidable. 

Proof. By Lemma [3 let us assume all eigenvalues have real part less than or 
equal to 0. Then, using Theorem [6l we may consider the system as being 
represented by 

m 

f{t)^Y.p=^y''- 

We may split this exponential polynomial in two (reordering as necessary) and 
write f{t) = fi{t) + /2(t), where 

fe 



fe 

E 



Pi,3 it) COS (A,t) + P2,, (t) sin(A,i)) (5) 



are those terms with real part and f2{t) is the summation of the remaining 
terms. Equation (O follows from Theorem [H] since fi{t) is real valued. Since 
{A, c, a;(0)) is reduced, fi{t) is not identically zero. Note that f2{t) tends to zero 
exponentially fast as t increases. 

For a polynomial P of degree n we may use Cauchy's bound on the maximum 
modulus of any polynomial root to determine that for any root z G C of P{x) = 
Qnx" + . . . + aix + ao we have that: 

max{|ao|,|ai|,...,|a„^-_i|} 

+ K\ ' 

as is easy to prove. Thus define T g M>o to be strictly greater than this 
maximum bound for any Pi j or P2,j in Equation ([5]) for 1 < j < m and thus 
for all t >T, the sign of Pij{t) and P2.j{t) for each 1 < j < m is fixed. 
For each 1 < j < k there exists tj^i, tj_2 > T such that 

■Pijfei)cos(Ajtja) +P2j(tj,i)sin(Ajtj, i) > 0, 

-Pijfe2)cos(Ajij-2) + P2j(ij,2)sin(Ajtj,2) < 0. 

Each Aj is distinct thus we have enough freedom in the choice of these times 
so that there exists ti , t2 > T such that 

Pij(ti) cos(Ajii) + P2j(ti) sin(Ajti) > , < • < , 
Pi,j{t2) cos(Aji2) + P2,j{t2) sin(Aji2) < ' ^-^ - 

We now see that /i(ti) is positive and /i(i2) is negative thus there exists a 
solution since /2 {t) decays exponentially fast and there exists an infinite number 
of solution times. □ 
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We now use Theorem [5] (Baker's theorem) to provide bounds on sums of 
exponentials polynomials. We first start with some lemmata. 

Lemma 12. Let wi and uji he different algebraic numbers, linearly indepen- 
dent over Q, and e*"^^ , e*"^^ be algebraic numbers on the unit circle. There exist 
effective constants C, N,T > such that at any time instant t > T , either 
1 - cos(cJit + (t>i) > C/t^ or 1 - cos{LU2t + ^2) > C/t^. 

Proof. We prove that there are C, N such that jwii + (pi — 2fc7r| > C/t^ for all 
integers k or \uJ2t + (j)2 — 2fc7r| > C/t^ for all integers k. 

Indeed, suppose that for some t, k, I and e > 0, both jojit + ^i — 2fc7r| < e and 
|<^2i+02— 2/7r| < e. Then \t+(l)i/uji—2kTr/uJi\ < e/|wi| and \t+(j)2/ 102— 21^/1^21 < 
e/\u!2\. By difference wc find I (/ii/tJi — (/)2/i^2 + 2Z7r/[j2 — 2fc7r/wi I < + )■ 

Let us introduce = (r^ + T^y^. Then\k—2TTi-l^2m+—i(j)i~^id>2\ < 
e. 

Observing that 2Tri,i(j)i,i(j)2 are logarithms of algebraic numbers, we apply 
Baker's theorem. Note that the height of ka, for any algebraic number a of 
degree d, is at most \k\'^ times the height of a, from the definition of height. 

Thus, e > max(Ci|fc|'*,C2|?|'',Co) °, for some Co,Ci,C2, Nq not depending 
on k, I. It is also clear that, given t, \ujit + (j>i — 2fc7r| is closest to zero for some 
k < C^t, and similarly for \uj2t + 4>2 — 2ln\. 

This proves that there exists C',Nq,T such that for every t > T, either 
\ujit + </)!- 2fc7r| > C't-^° for aU k, or \uj2t + 02 - 21t:\ > C't-^» for all I. 
Since 1 — cos(a + 2fc7r) > a^/S for a small enough and some k (from Taylor 
approximation), the claim follows. □ 

We say that a property T-eventually holds for a function g : K>o — *■ R if it 
holds for all time instants t > T. For instance, g is eventually positive if there 
is a threshold T such that g{t) > for all t > T. Clearly, if / is the solution of 
a linear differential equation, then it has finitely many zeros if and only if it is 
eventually positive or eventually negative. 

We say that gi is (T, r)- exponentially dominated by g2 if |5i(i)| < e^^*|(72(i)| 
for r > and aU t > T. 

Lemma 13. Let us consider T-eventually nonzero continuous functions gi, . . . , gk 
and the function f (t) = 5o(i)+X]j=i cos(ti;jt+(/)j), where uii,u!2 are linearly 
independent positive algebraic numbers, go is {T,r)- exponentially dominated by 
gi and g2, and (j>i,4>2 o-fe angles such that e*'^^,e*"^^ are algebraic. Then the 
following T' -eventually holds, for some T' : 

k k 

-J2\9j\ <f<J2\9j\- 
i=i i=i 

Moreover, such T' can be computed as a function of T,r,uJi,LU2, (t>i, 4>2- 
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// all LOj are linearly independent over the rationals, then for any e > 0, there 
exist arbitrarily large times t such that 

k 

I{t)>{l-e)Y,\g,{t)\ 
and arbitrarily large times t such that 

k 

f{t)<^{l-e)Y,\g,{t)\. 

i=i 

Proof. It is obvious that -|5o| - Ej=i Iffil < / < Iffol + ISj I- 

To prove that we can get rid of go, we exploit the fact that the cosines 
cos(tJit + (j)i) and cos(cl;2^ + 4>2) never take the value ±1 exactly at the same 
time, except possibly once; this is a consequence of the linear independence 
of wi and u)2. Due to Lemma [T^ one of the cosines is Ct~^ away from 1, 
for some C, N and all large enough times. Then for all large enough times 
t, either \go + giCos{LUit + (pi) + g2COs{uj2t + 4>2)\ < \gi\{l - Ct~^) + \g2\ or 
\go+gi cos{uJit + 0i) + g2Cos{uj2t + 02)| < l^il + 1521(1 - Ct^^). In any case, 
since go is exponentially dominated by both gi and g2, we have \go+gi cos{uJit + 
0i) + 92 cos{uj2t + 4>2)\ < \gi\ + I52I, for some T', computable as a function of 
T,r,C,N. Adding all the terms gj cos{u!jt + 4>j) proves the first claim of the 
theorem. 

We now prove the second claim. From Kronecker's theorem and the linear 
independence of frequencies, we have that the set F — {(cos(ijjji + 0j))i<j<fc|t > 
0} is dense in [—1, 1]*^. Hence, F will approach all the vertices of [—1, 1]'^' by less 
than any e > for arbitrarily large times. For those times such that for all j, 
cos{ujjt + <j)j) is close within e/2 to sign gj, and \go{t)/ gi\ < e/2, the first part 
of the second claim holds. The second part is similar. □ 

We now prove the main theorem of this section, which says that in some 
circumstances, one can reduce the search for a solution to an instance of Con- 
tinuous Skolem Problem to a finite time interval. Recall that an eigenvalue 
is nondefective if its algebraic and geometric multiplicities coincide. The fre- 
quency of an eigenvalue is the absolute value of the imaginary part. Recall that 
for a real matrix, complex eigenvalues come in conjugate pairs, determining one 
equal frequency. 

Theorem 14. Given an instance of Continuous Skolem Problem where 
all dominant eigenvalues are nondefective, at least four in number and such 
that the set of their distinct nonzero frequencies is linearly independent over the 
rationals. 
Then 

• The existence of infinitely many solutions is decidable; 
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• If there are finitely many solutions, then those solutions are in [0, T], where 
T is computable. 

Note that multiple dominant eigenvalues are allowed. 

Proof. As allowed by Lemma [3 we can suppose without loss of generality that 
the dominant eigenvalues are on the imaginary axis. 

Then we are looking for real zeros of a function f{t) = 70 + fi{t) + /2(i), 
where 70 is the contribution of the dominant zero eigenvalue (if any), 

fe 
fe 

= otj cos (Ajt) + I3j sin {\jt) 

collects the dominant terms corresponding to dominant complex eigenvalues 
9j = iXj and f2{t) is exponentially decreasing. By elementary trigonometric 
manipulations, /i can be converted into 

fc 

hit) ^ 'Y13 cos {Xjt + (I) j), 

for some such that exp(i0j) is algebraic. Hence /i is a linear combination 
of shifted cosines. 

Since there are at least four distinct dominant eigenvalues, /i contains at 
least two different frequencies. We apply Lemma [T51 with gi — ji,ga ~ f2 to 
obtain that the following eventually holds: 

fe fc 

Moreover, the same lemma tells us that for any e > there are arbitrarily 
large times t such that: 

fe 

-(l-e)^|7,|>/(t)-7o 
and arbitrarily large times t such that: 

fe 

/W-7o>(l-e)El7^l- 

As mentioned above, / has finitely many zeros if and only if / is even- 
tually positive or eventually negative. It results from the above that / < 
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(T-eventually, for some T) if and only if 70 + J2'j=i — 0- Moreover, when 
/ < (T-eventually, for some T), such a T can be be computed. A similar 
argument holds for / > 0. This proves the claim. □ 

Note that in discrete time, checking the existence of a zero in a finite time 
interval is a trivial task, while in continuous time we do not know how to decide 
the existence of a zero between time and T. 



5. NP-Hardness of Nonnegativity Problem 

We now prove the continuous version of Blondel-Portier's result 

Theorem 15. The nonnegativity problem for instances of CONTINUOUS Skolem 
Problem given by a skew- symmetric matrix is NP-hard and decidable in expo- 
nential time. In particular, the general nonnegativity problem is NP-hard. 

Proof. A skew symmetric matrix has only imaginary eigenvalues and Jordan 
blocks of size one. By Theorem |6] we must find nonnegative real zeros of a 
function of the form 

f{t) = X! cos(Ajt) -I- Pi sin(Ajt). 

i 

We can, in polynomial time, find a basis ^1 , . . . , over the rationals for the 
family Ai , . . . , , such that every Ai is an integral combination of ^1 , . . . , . For 
every we introduce two variables Xi = cos(^ii) and yi — sin(^it), which satisfy 
x'^ -\- yf = 1. Hence f(J,) is a polynomial P in Xi,yi (by elementary trigonom- 
etry). From Theorem [4] (Kronecker's theorem), the trajectory (^it, . . . ,£,kt) is 
dense in [0, 27r]'', from which {f{t)\t £ M} = {P{xi,yi, . . . ,Xm,ym)\xi,yi G K} 
follows. Hence, / is nonnegative if and only if P is, when taken over the set 
{xi,yi, . . . , Xm,ym\xj,yj S M and x'j -\- yj = 1 for 1 < j < to}. This problem is 
solvable in time exponential in to by Tarski's procedure (see for example Q)- 

Suppose we are given a polynomial P{xi, . . . ,Xk). We write Xi — cos(^it) 
for every 1 < i < k. Every monomial of P can therefore be written as a 
linear combination of cosines by elementary trigonometry. For instance, X1X2 = 
cos^iicos^2i = , and so on. In this way, the polynomial 

P{xi, . . . ,Xk) can be written as a function f{t) = cos(Ait), such that 

{/(t)|t e M} ~ {P{xi, . . . ,Xk)\xi G [—1,1]}. Hence / is nonnegative if and 
only if P is nonnegative on [—1,1]''. Since checking the nonnegativity of a 
polynomial on [—1,1]'^ is NP-hard (which is easily proved via an encoding of 
the 3-SAT problem, see, e.g., [lo|). then the nonnegativity problem for instances 
of Continuous Skolem Problem is also NP-hard. □ 

Note that physical linear systems that preserve energy can often be mod- 
elled by differential equations with a skew-symmetric matrix, because these are 
precisely, up to a change of variables, the systems for which the energy l/2x"^ a; 
(where x is the state) is constant along the trajectories, (see, e.g., [29|). This 
case is therefore particularly relevant. 
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6. Conclusion 

In studying this problem, we are not so much interested in exactly describing 
the solutions to the problem, as determining the existence of solutions. For 
example, if we have algebraic times ti,t2 G M>o with ti < t2 such that /(ti) 
and f{t2) have different signs then there exists t £ [ti, such that f{t) = by 
the intermediate value theorem. 

The main problem encountered in solving Continuous Skolem Problem 
however appears to be that f{t) can reach tangentially, i.e. we may have a 
solution f{t) — where there exists e > such that /(t) > for all r G [t — e,t+ 
e]. Since, by Lemma[51 the solution will, for non trivial cases, be transcendental, 
it is difficult to determine when such a situation arises. Indeed, given a real 
valued exponential polynomial, if we take its square then it is positive real valued 
and reaches zero tangentially if and only if the first exponential polynomial had 
a zero. 

We have therefore attempted to show several instances in which the problem 
is decidable but the general problem remains open. The equivalent problem of 
determining if an exponential polynomial has real zeros seems equally interest- 
ing. It is surprising that the problem is open even for a finite time interval. 
Solving Skolem's problem in the discrete case for finite time is obviously decid- 
able since we can simply compute the values in the interval. 

Open Problem 16. Is Bounded Continuous Skolem's Problem decid- 
able? I.e. Given a fixed T £ M>o, and an instance of Continuous Skolem 
Problem, f{t) = c'^exp(At)x(0), does there exist t <T such that f{t) =0? 

We also showed that the nonnegativity problem is NP-hard in the continuous 
case. It is not clear if a similar technique can be used to show that CONTIN- 
UOUS Skolem Problem is also NP-hard. In the discrete Skolem's problem it 
turns out that determining the nonnegativity and positivity of a linear recurrent 
sequence are equivalent in terms of complexity, however this is not clear in the 
continuous case. 

Open Problem 17. Are Continuous Skolem Problem and the continuous 
nonnegativity problem computationally equivalent? 
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